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Résumé : 



La microstructure des aciers dual-phases et leur comportement mécanique sont fortement sensibles 
aux variations de procédé (traitements thermiques). Un contrôle en ligne par méthode magnétique est 
^/^ • envisagé, la mesure s 'appliquant à un matériau soumis à un état de contrainte. Le dual-phase est 
Q.^ . un milieu biphasé (ferrite/martensite), où chacune des phases est considérée comme étant une sphère 
immergée dans un milieu homogène équivalent. La modélisation de chaque phase s 'appuie sur un modèle 
magnéto-mécanique couplé. Il s'agit d'un modèle monocristallin explicite représentatif du poly cristal 
pi |! isotrope de la phase considérée. La mise en place de règles de localisation permet la simulation du 

' milieu biphasé. Expériences et modèle sont comparés. 
c/3 ! 

^ ■ Abstract : 
O \ 

^ , The microstructure and mechanical behavior of dual-phase steels are highly sensitive to the variation 
Ô ' of the process (heat treatments). Online control by magnetic method is relevant. A measurement under 
'S . applied stress must be considered. The dual-phase is a two-phase médium (ferrite / martensite). Each 
phase can be considered as a sphère embedded in a homogeneous équivalent médium. The model used 
for each phase is based on a magneto-mechanical coupled model. This is an explicit single crystalline 
model représentative of the behavior of the corresponding phase. Localization rules allow the simulation 
of the two-phases médium. Experiments and modeling are compared. 

• Mots clefs : comportement magnéto-mécanique, microstructure biphasée, localisation. 

1 Introduction 

■ Les dernières années ont vues un intérêt grandissant des industries automobiles pour l'utilisation 
d'aciers à haute performance tels que les aciers dual-phases (DP). La production de ces aciers implique 
plusieurs procédés de fabrication : métallurgie primaire et métallurgie secondaire, coulée, laminage (à 
chaud/froid) et traitements thermiques. Ces procédés conduisent à une microstructure biphasée prin- 
cipalement composée d'îlots de martensite dure dispersés dans une matrice ferritique ductile (figure [T]) 



^ ■ en proportion variable selon l'histoire thermo-mécanique du matériau. Leur microstructure (fraction 



et composition des phases) étant fortement sensible aux variations du procédé de fabrication (traite- 
ment thermique, laminage). Les industriels cherchent à mettre en place une méthode de contrôle non 
destructive en ligne. La méthode employée exploite les propriétés magnétiques de ces matériaux. La 
mise en œuvre de ce procédé de contrôle rentre dans le cadre du projet ANR DPS-MMOD [T]. 

Le travail présenté dans cette communication a été réalisé dans le cadre de ce projet : il s'agit d'une 
contribution à la modélisation du comportement magnéto-mécanique des aciers DP. La méthode de 
modélisation retenue s'appuie sur le modèle multidomaine, développé au LMT-Cachan. Ce modèle 
est basé sur la partition en domaines magnétiques des matériaux magnétiques. Il s'agit d'un modèle 
monocristallin. On montre sous certaines hypothèses qu'il existe une direction de sollicitation pour 
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Figure 1 - Microstructure d'un acier dual-phase : martensite en blanc, ferrite en sombre. 



laquelle le comportement obtenu est représentatif du polycristal isotrope du matériau concerné. Il 
s'agit, dans le cadre de l'application recherchée, de prendre en compte : 

- l'orientation du chargement magnétique dans le repère du monocristal représentatif. 

- un chargement mécanique statique imposé et ses effets sur la mesure. 

- la nature biphasée de la microstructure de l'acier DP, ce qui suppose de mettre en place des règles 
de localisation mécanique et magnétique adaptées. 

2 Modélisation multidomaine |j2|, [3] 

Les deux phases en présence (ferrite / / martensite m) sont ferromagnétiques : elles s'aimantent et 
se déforment en présence d'un champ magnétique. Les deux phases sont modélisées séparément à 
l'aide d'un modèle de comportement proposé récemment : le modèle multidomaine [2]. Ce modèle 
prend son origine dans la configuration en domaines magnétiques des matériaux magnétiques. On 
considère un monocristal de fer de symétrie cubique divisé en 6 familles de domaines notés a ( axes 
< 100 > du cristal) correspondant aux 6 directions de facile aimantation (figure [Jh-) ■ On note 7" = 
(7° 72 72 1^ vecteur directeur de chacun des domaines ([T]). Il est colinéaire au vecteur aimantation 
M" = où Mg désigne l'aimantation à saturation du matériau. Le vecteur aimantation est 

initialement confondu avec les axes du cristal. Sous l'effet d'un chargement magnétique ou mécanique, 
les fractions de domaines ainsi que les directions d'aimantation changent. 




Figure 2 - (a) représentation simplifiée du monocristal ; (b-c) triangle standard et paramétrage du 
chargement. 

Si on utilise des hypothèses de déformation et de champ homogène sur le cristal [1], l'énergie libre 
totale par domaine est la somme de trois contributions énergétiques : 
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- L'énergie de champ : = —fioH .M" . 

- L'énergie magnéto-mécanique : = —cr : ej^. 

- L'énergie magnétocristalline : = Ki (717I + 7273 + 7173) 

est la perméabilité du vide (47r.l0~^ H/m), H et a sont le champ magnétique et le tenseur des 
contraintes appliqués au cristal, e° est le tenseur de magnétostriction du domaine considéré ([2]), K\ 
est la constante d'anisotropie principale du matériau. 



Aioo(7i - \) Aiii(7i72) Aiii(7i73) 



' Aioo(7|-è) Ani(7273) 1 (2) 



Il 2 



Tir 



sym Aioo(73 - \] 

( cos (/>c sin 9c sin Oc sin 4>c cos ) (3) 



Compte tenu de la symétrie cubique du cristal, une direction de chargement ( champ magnétique ou 
traction uniaxiale) trouve son équivalent dans le triangle standard défini par les axes [100], [110] et 
[111] du cristal (figure [2)o). On considère une direction de chargement fie définie par les deux angles 
sphériques (j)c et Oc ([3]). Les directions de chacun des domaines minimisent par définition l'énergie libre 
totale. Compte tenu de la restriction au triangle standard, il est possible de réaliser une minimisation 
analytique. On aboutit par exemple aux expressions suivantes pour le domaine (1) de la figure [2^ : 

1,^ , _ HoMgHcjyc + arctan (|Aiiio- sin(2(jic)) 
^ ' ^' ~ HoMsH + 2Ki + 3A100CT cos(20c) 

(4) 

1 = ^o^^^d - + arctan (f Ama sin(2(f - Oç)) 

{ ,(^,nc) - ^ HoMsH + 2Ki + 3A100CT cos(2(f - Oc)) 

avec H = H.fîc et a = a.fîc les vecteurs champ magnétique et contrainte. Les fractions de chaque 
domaine /" sont ensuite obtenues en utilisant une relation explicite inspirée d'une fonction statistique 
de Boltzmann U] : 

^ exp{-A,.W-) 

Ea expi-As.W^) ^ ' 

As est un paramètre proportionnel à la susceptibilité initiale xo du matériau tel que Ag = ^J^j2 ■ Le 
comportement étant considéré homogène dans le monocristal, les règles d'homogénéisation s'appliquent 
pour le calcul des valeurs moyennes de l'aimantation et la déformation : 



On aboutit à un modèle capable de rendre compte du comportement magnéto-mécanique d'un mono- 
cristal. Or, puisque toutes les directions de chargement possibles peuvent être restreintes au triangle 
standard, le comportement d'un polycristal isotrope est nécessairement donné par un chargement le 
long d'une direction appartenant à ce triangle. La non linéarité des phénomènes en présence implique 
que cette direction n'est pas une direction moyenne, et est susceptible de changer d'orientation avec la 
contrainte ou le niveau du champ magnétique imposé. On fait l'hypothèse que ce changement est assez 
petit pour être négligé [3]. La direction moyenne théorique est obtenue pour (pc = 39° et Oc = 78° [5] . 

3 ]V[odélisation d'un acier DP, localisation des champs 

La présence de plusieurs phases de nature différente crée une perturbation locale du champ (magnétique 
ou mécanique) imposé à chacune des phases. Le champ local ne vaut ainsi généralement pas le champ 
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moyen. On parle de localisation du champ macroscopique considéré. On considère le milieu biphasé 
{f,m). On note / une inclusion, jouée alternativement par la phase / ou la phase m. 

Le champ magnétique local Hj vu par la phase / est une fonction complexe du champ macroscopique 
H et des propriétés du milieu moyen. Dans le cas d'un problème d'inclusion sphérique [3], et en utilisant 
une hypothèse de champ homogène par phase et de comportement linéaire, la relation de localisation 
prend la forme : 

Xo et M représentent la susceptibilité et l'aimantation du milieu moyen, Mj étant l'aimantation de 
l'inclusion. L'extension au comportement non linéaire suppose d'utiliser une susceptibilité sécante pour 
la définition de xoi soit xo = ll-^ll/ll-f^ll- Cette approche est appliquée à la microstructure biphasée 
{f,m). On obtient : 

Hf = H + ^[M-Mf) ; H^ = H + ^[M-Mm) (8) 

avec : 

H = ffHf + fmHm ; M = ffMf + fmMm (9) 

Les champs Hf et Hm sont introduits en entrée de deux calculs multidomaines. Ces calculs fournissent 
une aimantation par phase, puis une moyenne permettant une nouvelle évaluation des champs locaux. 
Ce processus est itéré jusqu'à stabilité des grandeurs magnétiques (illustré figure H^). 

La solution du problème d'inclusion d'Eshelby constitue la base de modélisation du comportement 
des milieux hétérogènes en mécanique. La déformation de magnétostriction Ej est la déformation que 
subirait l'inclusion considérée en l'absence de la résistance exercée par la matrice. Aux contraintes 
d'incompatibilité s'ajoute une contrainte associée au chargement extérieur a. Si de plus on décompose 
les déformations en une somme de déformation élastique et de déformation de magnétostriction et qu'on 
considère les constantes d'élasticité homogènes (même constantes élastiques pour les deux phases), on 
aboutit à la formulation suivante de la contrainte crj dans l'inclusion / : 

aj = a + C{I- S^) : (e^ - e^) (10) 

C est le tenseur d'élasticité du milieu ; E>^ est le tenseur d'Eshelby, ne dépendant que des paramètres 
matériau de la matrice et de la forme de l'inclusion, e'^ est le tenseur de magnétostriction moyen. Cette 
approche est appliquée à la microstructure biphasée (/,m) en considérant des modules et des formes 
identiques. On obtient : 

crj = cr + C(I - S^) : (e^^-e^) cr^ = cr + C(I - S^) : (e^ - e^,) (11) 

avec : 

a- = fjaf + fmCTm ; = fje'^ + fm€m (12) 

Les champs ctj et a m sont introduits en entrée de deux calculs multidomaines au même titre que 
les champs magnétiques. Ces calculs fournissent une déformation et une aimantation par phase, puis 
des moyennes permettant une nouvelle évaluation des champs mécaniques et magnétiques locaux. Ce 
processus est itéré jusqu'à stabilité des grandeurs. Restant dans un cadre mécanique linéaire, la mise 
en place numérique de cette deuxième localisation ne pose pas de problème majeur. 
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4 Résultats expérimentaux et modélisation 

Une première étape a consisté à identifier les paramètres du modèle multidomaine pour les deux phases. 
Nous disposons d'échantillons de fer pur et d'un acier XClOO trempé à l'eau dont la structure est 100% 
martensitique. La figure [S] regroupe les résultats expérimentaux obtenus pour les deux matériaux (com- 
portements magnétique et magnétostrictif). On se reportera à [5] pour les procédures expérimentales. 
Les simulations à l'aide du modèle multidomaine monophasé sont présentées sur les mêmes figures. Le 
tableau [T] regroupe les paramètres utilisés. Ceux du fer sur sont issus de la littérature. Les paramètres 
de la martensite sont mal documentés. Nous avons procédé à une optimisation pour leur évaluation. 
Le coude de saturation est relativement mal décrit par le modèle pour les deux matériaux. C'est un 
défaut associé à l'hypothèse de champ homogène intrinsèque au modèle multidomaine. 
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Figure 3 - Comportements magnétique (a,c) et magnétostrictif (b,d) du fer pur (a,b) et de la mar- 
tensite pure (c,d) - comparaisons modèle/expérience. 
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Aiii 
-21.10" 
3.10-6 



Ki{J.m~'^) 
4,8.10^ 
10.10^ 



1,71.10*^ 
1,05.10*^ 



3,5.10-3 
4.10-^ 



Table 1 - Paramètres utilisés pour la modélisation. 



Une fois les paramètres du comportement de la ferrite et de la martensite identifiés, il est possible de 
simuler le milieu biphasé. On utilise une fraction de 42% de martensite (fraction identifiée à partir 
des micrographies). La figure H] permet d'apprécier l'effet de la procédure de localisation en champ (le 
calcul se fait en contrainte homogène) . La figure illustre la convergence des champs locaux après 
quelques itérations : la phase magnétiquement dure (martensite) voit un champ plus élevé que le champ 
moyen ; le champ perçu par la ferrite est plus faible que le champ moyen. La figure [Hd montre l'effet de 
cette localisation sur le comportement magnétique. On constate que l'aimantation moyenne prévue est 
plus faible avec localisation. La figure [5] nous montre le résultat des simulations comparé aux résultats 
expérimentaux. Comme on peut le voir, le comportement magnétique est bien décrit par le modèle, 
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passé 600A.m~ , on voit que la procédure de localisation n'a plus qu'un effet marginal. On constate que 
la déformation de magnétostriction est surestimée à aimantation élevée. La procédure de localisation 
en contrainte devrait permettre d'améliorer ce résultat. 




Figure 4 - Illustration de la localisation (/m=42%) - (a) convergence des champs locaux pour 
H = 1,25.10^ A.m~^ -(b) Comportement magnétique : expérience et simulation. 




H (A/m) M (A/m) j^^gS 

Figure 5 - Comparaison modèle/expérience pour le comportement magnétique (a) et magnétostrictif 
(b) du dual-pliase (/^=42%) 



5 Conclusion 

Le modèle proposé permet un estimation rapide du comportement magnéto-mécanique d'une micro- 
structure biphasée. Cette estimation requiert de connaître le comportement de chacune des phases. 
On constate une surestimation du comportement magnétostrictif, imputable à la non prise en compte 
de la localisation en contrainte ou à une modélisation trop simpliste (magnétostriction isotrope) de la 
martensite. La technique développée servira à faire dialoguer des méso- modèles micromagnétiques [T]. 
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